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^ ' Double-tagged interactions of photons with virtualities Q^ between 10 GeV^ 

O ■ and 200 GeV^ are studied with the data collected by DELPHI at LEPll from 

^ I 1998 to 2000, corresponding to an integrated luminosity of 550 pb~^. The 
7*7* -^ yU^/^^ data agree with QED predictions. The cross-section of the 
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reaction 7*7* — > hadrons is measured and compared to the LO and NLO 
BFKL calculations. 
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1 Introduction 

This paper presents the study of double-tagged two-photon interactions e^e^ -^ 
e"''e~7*7* -^ e^e~X (where X is either a muon pair or hadrons) with the DELPHI 
detector [1] at the CERN LEPII colhder. Both scattered electrons^ are detected by the 
Small angle Tile Calorimeter (STIC). Compared to the untagged or single-tagged modes, 
with both or one of the electrons escaping detection, this mode of gamma-gamma collision 
has the advantage that the kinematics of the interaction is well defined by the measure- 
ment of the energies and scattering angles of the tagged particles. The production of 
muon pairs is described by QED. Similarly, multihadron production is expected to be 
described by QPM, but only in a first approximation. If the virtualities of the photons 
are large enough, it is predicted that there should be a large contribution from processes 
with (multi)gluon exchange between qq dipoles [2], which is described by the BFKL equa- 
tion [3]. Two-photon interactions are therefore a suitable process to investigate BFKL 
dynamics. Figure 1 shows the main diagrams relevant to the analysis. 

The kinematics of the process is illustrated in Figure 2. We use the following notations: 
Pi (i=l,2) are the four-momenta of the beam electrons, y^ is the e^e~ centre-of-mass 
energy, Ebeam is the beam energy; the four-momenta of the scattered electrons, their polar 
angles and their energies are p^, 9i and Ei respectively. 

The variables relevant to this study are the virtualities of the photons, Qf, the invariant 
mass of the two photons W^*^* and a dimensionless variable Y: 

• Ql = -qf = -{P^-P^? = 4E,E,,,^sm\e,/2); 

• W^,^, = -(gi + q^f ~ syxy2 with t/i = 1 - {Ei/E^^am) cos^iOi/2); 
.F = ln(l^^%./v/Qfgf). 

The Y variable is used to compare the multihadron data with the BFKL predicted 
cross-section with the conditions W^,^, ^ Qf and | \n{Ql/Ql) |< 1, where the second 
condition is needed to select virtualities of the photons of the same order. 

The analysis is divided into two parts: the study of the production of muon pairs aims 
at comparing the data with the well-known QED model and at tuning the experimental 
cuts, while the multihadron production is used to measure the cross-section a^.^. and 
to compare it with the BFKL predictions. The models used for each part and the back- 
ground estimations are described separately. 



2 Detector and data sample 

A detailed description of the DELPHI detector and of its performance is presented 
in ref. [1]: here only the components relevant to the present analysis will be briefly 
mentioned. 

The scattered electrons are detected in the luminosity monitor STIC, which covers 
the region from 29 mrad to 185 mrad in the polar angle 6 ^, with R(f) segmentation of 3 
cm x22.5° [4]. Given the energy and angular resolution of the STIC calorimeter, the Q^ 

'^Throughout this paper, electron stands both for electron and positron. Asterisk over 7 symbol explicitly indicates that 
the photon is highly virtual. 

■^The origin of the DELPHI reference system was at the centre of the detector. It coincides with the ideal interaction 
point. The z-axis was parallel to the e~ beam, the x-axis pointed horizontally to the centre of the LEP ring and the y-axis 
was vertically upward. The co-ordinates R, (j>,z formed a cylindrical coordinate system and was the polar angle with 
respect to the z-axis. 



resolution varies between 1 GeV^ and 2.5 GeV^ in the Q^ domain of the present analysis 
(gf between 10 GeV^ and 200 GeV^). 

Charged particles are detected in the barrel tracking system comprising the Silicon 
Tracker (ST), the Inner Detector (ID), the Time Projection Chamber (TPC), and the 
Outer Detector. In the endcap regions, they are detected by the ST, by the TPC down 
to 20° in polar angle, by the ID down to 15° and by the Forward Chambers A and B. 
All detectors are located inside a superconducting solenoid providing a uniform magnetic 
field of 1.23 T parallel to the axis of the colliding e^e~ beams. The combined momentum 
resolution provided by the tracking system is a few per-mille in the momentum range of 
this study. 

Muon tagging is performed with the Barrel and the Forward muon drift chambers, 
and with the Surround Muon Chambers based on limited streamer tubes, which cover 
the gap between the previous two. 

The study is done with the DELPHI data taken during 1998-2000 at e^e~ centre-of- 
mass energies from 189 GeV to 209 GeV, corresponding to an integrated luminosity of 
550 pb~^, with the subdetectors relevant for the analysis all fully operational. 

Simulated events for the physics processes and backgrounds are generated at the dif- 
ferent centre-of-mass energies and passed through the full DELPHI simulation and re- 
construction chain. 

3 Study of e+e~ — ^ e+e~ + /x+/x~ interactions 

3.1 Data analysis 

The following criteria are used to select 7*7* -^ fi^fi^ events: 

• There are two clusters with energy deposition Ei greater than 30 GeV, one in each 
arm of the STIC ^, and the polar angle 6i (or (180°-^^) if the angle is more than 90°) 
exceeds 2.2° for each cluster; 

• Qf is between 10 GeV^ and 200 GeV^ for both tagged particles; 

• The acollinearity of the scattered electrons is above 0.2 degrees. This criterion 
removes a superposition between Bhabha events and untagged 77 -^ fi~^fi~ events; 

• Each event contains two charged particles with zero net charge and invariant mass 
between 2 GeV/c^ and 50 GeV/c^. Particles are considered if their momentum is 
greater than 400 MeV/c, their polar angle is within the interval 20° - 160° and their 
impact parameters are smaller than 4 cm in R(f) and 10 cm in z; 

• At least one of the charged particles is identified as a 'standard' or 'tight' muon by 
the DELPHI tagging algorithm [1]. 

The number of selected e+e^ -^ e+e^ + fi^fi^ events is 226. Double-tagged events 
are triggered either by the STIC trigger component or by a single charged particle track 
component [5]. The trigger efficiency has been calculated using the redundancy of the 
trigger together with independent calculations based on the parameterization of the single 
track efficiency [5], and is found to be larger than 99%. 

The BDKRC event generator [6], including the full set of QED diagrams and all 
fermion masses, is used for the Monte Carlo simulation of 7*7* -^ ^^ ^~ events. The 
processes corresponding to diagrams other than the diagram of two-photon interactions 

■^The choice of the cutoflf for the ininimum energy of the tagged particles has to be made carefully. It should be as small 
as possible since low electron energy corresponds to large values of VKy*-,* and Y (which are important in the multihadron 
case, see below). At the same time the lowering of this cut leads to an increase of the off- momentum background and thus 
decreases the accuracy of the measurement. 



(multiperipheral) are found to give a contribution of about 2 percent. The number of 
multiperipheral events is expected to be 194. 

The following sources of background are considered: 

• The coincidence of a gamma-gamma event with an off-momentum electron. The 
probability of such coincidences, averaged over the data from different years, is 
calculated with 77 -^ fi'^fi~ events to be (0.0016 ±0.0002). Using this value, the 
background from a superimposition of two off-momentum electrons with an untagged 
77 — s> fj,~^fj,~ event turns out to be negligible; 

• The coincidence of one off-momentum electron with a 7*7 — > fi^fi~ single-tagged 
event, i.e. when one scattered electron is detected in the STIC while the other one 
is an off-momentum electron. Usually this background is evaluated by convoluting 
the 7*7 -^ fi'^fi~ Monte Carlo simulation with the spectrum of off-momentum 
electrons. This requires the appropriate description of the single-tagged data by 
the simulation and thus it is model-dependent. The following approach avoids the 
problem and calculates the background directly from the data. The cuts as listed 
above are applied with one difference: events with two electrons in the same STIC 
arm and none in the other are selected. They include one off-momentum electron. 
Then one of the electrons is rotated to the "empty" STIC arm, i.e. its pz component 
is inverted. The background coming from the coincidence of one off-momentum with 
a ry*ry _> ^+^~ sluglc-tagged event is thus estimated to be (15 ± 4) events; 

• The background from 7*7* -^ t^t~ events is estimated as (23 =b 2) events by using 
the TWOGAM event generator [7]. 

The overall background is thus estimated as (38 ± 4) events. 

Figure 3 shows the distributions of the tagged particle energies*^ normalized to the 
beam energy, their polar angles (two entries per event for both histograms), the invariant 
mass of the muon pair, W^^, calculated from the muon 4- momenta, and the distribution 
of the normalized longitudinal momentum balance defined as 

NLMB =1 PzMgl + Pz,ta92 + Pz,X I /Ebeam, 

where Pz^x is the z— component of the momentum of the system produced in the 7*7* 
collision. This variable has to be peaked at zero for well-reconstructed events and it is 
sensitive to the final state radiation. 

The analysis of 7*7* — > fi~^fi~ events permits the quality of the reconstruction to be 
examined. The kinematics of the gamma-gamma system is completely determined by the 
measurements of the muons and of the tagged electrons in this exclusive channel, and 
some quantities can be calculated either from the tagged particles or from the muons. 
Figure 4a (b) shows the difference between the gamma-gamma invariant mass {Y vari- 
able) calculated from the muon 4- momenta W^^{Y^^) and that reconstructed from the 
tagged particles' measurements W^^^iY^^). The asymmetry of the distributions is due 
to radiative corrections, as it has been verified by Monte Carlo simulation. These com- 
parisons show that the use of tagged particles is a good approximation to calculate the 
kinematic variables of the gamma-gamma system. 

3.2 Results 

The selected data sample is used for the measurement of the cross-section (o"ee ) of 
the reaction e"'"e~ -^ e~^e~ + /i"'"/i~ . The corrections for the detector acceptance and 
efficiency are done with the BDKRC simulated events. The statistical uncertainty in the 

''If the measured value of the electron energy is greater than -Ei,eami it is changed to (£^6eam"0-5 GeV) to be able to 
calculate the 7*7* invariant mass. 



MC simulation is included in the systematic error. Additional systematic uncertainties 
are evaluated by varying the selection criteria on the tagging particles. Systematic un- 
certainty coming from the muon identification procedure is negligible (second item of [1], 
p. 96 and references therein). 

The measured total o"ee cross-section is (1.38 ± 0.12{stat) ± 0.06{syst)) pb for virtu- 
alities of the interacting photons, Qf, between 10 GeV^ and 200 GeV^ and for invariant 
mass W^^ between 2 GeV/c^ and 50 GeV/c^. The QED expectation, including radiative 
corrections, is (1.36 ±0.01) pb. The cross-section calculated without radiative corrections 
is about 8% lower. 

The (Tee cross-section can be expressed via the fiux of photons with different polar- 
ization and the corresponding partial cross-sections of the 7*7* -^ ^~^fi~ interaction. 
This extraction procedure is described in the next section for multihadron production. 
Here only the result for the differential 7*7* — > fi'^fi~ cross-section as a function of Y 
is shown in figure 5. There is good agreement between the measurements and the QED 
predictions. 

4 Study of e+e~ — ^ e+e~ + hadrons interactions 

4.1 Data analysis 

The selection criteria for tagging electrons and for the charged particles are the same 
as described in the previous section. The sample of 7*7* -^ hadrons events is then 
selected by the following criterion: 

• Each event contains at least 3 charged particles with the invariant mass calculated 
from the particles' 4-momenta, Whad, larger than 2 GeV/c^; 

The following additional cuts are applied to suppress background events: 

• If the energy of one cluster in STIC, normalized to the beam energy, is larger than 
0.85 then the energy of another cluster has to be below 0.5. This cut is intended to 
suppress the contamination coming from e"'"e~ -^ hadrons events; 

• The thrust value of the charged particles, calculated in their centre-of-mass system, 
is less than 0.98 for the events with charged multiplicity below 5. The cut removes 
most of the 7*7* -^ t^t~ events. 

After these requirements, 434 events have been selected. Again the trigger efficiency [5] 
can be estimated from the redundancy of the trigger and from a parameterization of the 
single track efficiency, and is larger than 99%. 

The event generators used to simulate the 7*7* events and the background processes 
are listed below as well as the respective expected contributions: 

• TWOGAM (version 2.02) [7] and PYTHIA (version 6.205) [8] event generators (both 
include radiative corrections) are used to simulate 7*7* interactions. The expecta- 
tions are (331 ± 8) and (330 ± 8) events, respectively. The Monte Carlo generators 
include the quark-parton model (QPM) part and also the leading-order predictions 
for the resolved photon contribution; 

• The background coming from the process e^e~ -^ hadrons is simulated with the 
KK2f generator (version 4.14) [9] and its contribution is estimated to be (27 ± 3) 
events; 

• The background of r pairs produced in e^e~ annihilation is found to be negligible; 



• The contamination of r pairs produced in the two-photon interactions is evaluated 
as (26 ± 3) events by using the TWOGAM program; 

• The coincidence of an off-momentum electron with a 7*7 -^ hadrons single-tagged 
event is evaluated as (5 ± 2) events by using the same approach as described in the 
previous section. 

The data distributions for the photon virtualities, Qf (two entries per event), the 
invariant mass of the hadron system calculated with the charged particles' 4-momenta, 
Whad, the charged particles' multiplicity and the Y variable calculated with W^*^* are 
compared with the Monte Carlo simulation in figure 6. The data are represented with 
error bars. The solid and dashed histograms correspond to the sum of 7*7* simulated 
events obtained with the PYTHIA and TWOGAM generators, respectively, and various 
background sources. The hatched histograms show the estimated background contamina- 
tion. Both 7*7* models agree reasonably well with the data. The excess of the data over 
Monte Carlo (for low Q^ , large W) already indicates that the QPM term is insufficient. 

The calculations of the detector acceptance and efficiency have been done for both 
7*7* models and are shown in figure 7. The detection efficiencies express slightly different 
behaviour - the TWOGAM values are larger than the PYTHIA ones for high values of 
the Y variable, while for small Y values the behaviour is the opposite. The decrease of 
the efficiencies for Y above 4 is due to the selection criteria. 

4.2 Results 

The background subtracted data are corrected for detector effects using the two mod- 
els, and the measured differential cross-sections dcjee/dY are shown in figure 8 together 
with the average expectation of the two event generators used. The uncertainty due to 
the migration of events caused by the finite Q^ resolution (the relative uncertainty is 
around 0.08) is found to be small in comparison with the statistical uncertainty of the 
measurement. Note that, irrespective of the model, the data indicate a somewhat larger 
cross-section compared to expectations for high values of Y corresponding to large invari- 
ant masses of the 77 system. The total cross-section o"ee of the e+e~ -^ e~^e~ + hadrons 
interactions, within the phase space limited by the criteria Qj between 10 GeV^ and 200 
GeV^, and Whad above 2 GeV/c^, is measured to be (2.09 ± 0.17) pb using the correc- 
tions for detector effects based on TWOGAM and (1.86 ± 0.14) pb for the corrections 
based on PYTHIA. The statistical and systematic (see later) uncertainties are added in 
quadrature. The expectation of the quark-parton model is (1.81 ± 0.02) pb as obtained 
with TWOGAM. 

The 7*7* -^ hadrons interactions are expected to be sensitive to multiple gluon 
exchange (fig.l). The multigluon ladder is described by the BFKL equation [3], which 
predicts a growth of the cross-section at large Y. Note that the BFKL calculations are 
valid provided W^,^* ^ Qf (the variable Y should be larger than 2) and | \n{Ql/Ql) \< 1 
(to maintain the photon virtualities approximately equal). The application of this latter 
condition has the effect of reducing the data sample by about 37%. It has to be mentioned 
however that the migration of events around the chosen cut at unity does not introduce 
an appreciable systematic uncertainty. According to the Monte Carlo simulation, around 
3% of the selected sample had the true value of the logarithm above unity and pass the 
cut due to the Q^ resolution. Approximately the same percentage migrates inversely. 

The experimental conditions of the present study (Qf S> rn^ and the symmetry re- 
quirement for tagged particle detection) permit the relation between o"ee and a^^ , 
which initially reads [10] (the interference terms are omitted): 



(Tf. 



l^i,j=T,L ^ij^ij-: 



to be simplified to a relation involving an effective cross-section a^*^* , 

(Tee = LTT(Jrr "^^^^ ^I'l* = ^TT + 2e(TLT + (^^CTLL, 

where L^t is the flux of the transversely polarized photons calculable in QED, e is around 
0.94, aiT - 0.2aTT and ctll ~ O.OScttt [H]- The TWOGAM event generator including 
QED radiative corrections has been used to calculate Ltt'- it uses the decomposition of 
the cross-section for different photon helicities [10]. The limits on Qj, \ ln{Ql/Q2) \ and 
Whad are the same as described above. The differential cross-sections dcjee/dY , both for 
data and MC, and da^*^*/dY are all presented in Table 1. The selection efficiency is 
calculated using the mean of the results using TWOGAM and PYTHIA. The difference 
between the results using these two generators is used to calculate the systematic error due 
to modeling and is in included in the quoted systematic uncertainties. The flux binned 
over Y is presented as well. The absolute uncertainties on the dLxr/dY calculations are 
of the order of 0.002 x 10"^ 
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0.17±0.04(stat)±0.02(syst) 


0.26±0.03 


0.113 


0.89±0.23(stat)±0.21(syst) 


0-1 


0.42±0.06(stat)±0.08(syst) 


0.48±0.01 


0.090 


2.50±0.49(stat)±0.58(syst) 


1-2 


0.41±0.04(stat)±0.01(syst) 


0.48±0.01 


0.082 


3.56±0.42(stat)±0.08(syst) 


2-3 


0.30±0.03(stat)±0.02(syst) 


0.32±0.01 


0.070 


3.00±0.33(stat)±0.19(syst) 


3-4 


0.25±0.02(stat)±0.01(syst) 


0.15±0.01 


0.054 


2.83±0.32(stat)±0.02(syst) 


4-5 


0.08±0.01(stat)±0.01(syst) 


0.06±0.01 


0.034 


1.47±0.29(stat)±0.07(syst) 


5-6 


0.09±0.03(stat)±0.02(syst) 


0.02±0.01 


0.015 


3.51±1.33(stat)±0.60(syst) 



Table 1: The measured and expected differential cross-sections {daee/dY)data of the 
reaction e~^e~ -^ e^e~ + hadrons , the calculated photon flux Ltt including the radiative 
corrections [7], and the measured cross-section da^*^* jdY of the process 7*7* — * hadrons 
are shown as a function of the variable Y . 

The measured differential cross-section da^*^* / dY is shown in figure 9. The sys- 
tematic uncertainties are dominated by the difference between the results obtained with 
TWOGAM and PYTHIA Monte Carlo generators. The other systematic uncertainties, 
calculated by varying the selection criteria, the Q^ domain, etc. represent between 17% 
and 29% of the statistical uncertainties. The predictions of the QPM and of BFKL cal- 
culations, both in LO [12] and NLO [13], are also shown in figure 9. The two curves 
for the BFKL calculations correspond to the Regge scale parameter sq, which defines 
the start of the asymptotic regime, equal to Q^ or AQ"^. The LO calculations are more 
affected by the uncertainty coming from the choice of the scale parameter than the NLO 
ones. Note that the BFKL calculations are weighted over a number of Q^ bins and that 
therefore the running of a^ is also included. The data lie in any case much lower than 
the BFKL cross-sections calculated in leading-order. On the other hand, the data are 
closer to NLO predictions, since the expected growth of the gluon exchange contribution 
(BFKL) is much weaker and appears mainly for Y values larger than 4. Below this value 
the cross-section is dominated by the decrease of the QPM contribution. Unfortunately 
the LEP energy and the present statistics are not sufficient to study in detail the region 
at large F, where BFKL is expected to dominate. 



5 Conclusions 

Double-tagged 7*7* interactions have been studied with the DELPHI data taken 
at e+e~ centre-of-mass energies from 189 GeV to 209 GeV and corresponding to an 
integrated luminosity of 550 pb~^. For virtualities, Q^, of both photons between 10 
GeV^ and 200 GeV^ and final state invariant mass W above 2 GeV/c^, the cross-section 
of the process e+e^ -^ e+e^ + fi^fi~ is measured to be (1.38±0.12(stat) ± 0.06{syst)) 
pb, to be compared with the expectation of (1.36±0.01) pb for the QED calculations 
including radiative corrections to the photon fiux. The cross-section o"ee of the e+e^ -^ 
e^e~ + hadrons interactions is measured to be (2.09 ± 0.09(stat) ± 0.19{syst)) pb with 
the corrections for detector effects based on the TWOGAM event generator [7]. The 
differential cross-section da^*-y*/dY of the 7*7* — ^ hadrons interactions is measured 
and is compared with the predictions based on LO and NLO BFKL calculations. The 
leading order calculations clearly disagree with the data while the next-to-leading order 
predictions are found to be more consistent with the data, although the LEP energy is 
not sufficient to see a sizable effect due to the BFKL type contribution. The DELPHI 
data are in agreement with the results of the other LEP experiments [14]. 
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Figure 1: Main diagrams corresponding to the 7*7* -^ hadrons process. 
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Figure 2: The kinematics of 7*7* interactions. 
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Figure 3: (a) Normalized tagged particle energy Ei/Ebeam, (b) tagged particle polar 
angle 6i (two entries per event for both), (c) invariant mass of the muon pair W^^^ and 
(d) NLMB variable. The data are shown with error bars; the histograms are the sum 
of the BDKRC simulated events and of the estimated background. The hatched areas 
represent the background contamination. 
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Figure 4: (a) Difference between the invariant mass of the muon pair, W^^, and the 
invariant mass calculated from the tagged particles, (b) The same for the Y variable. 
The data are shown with error bars; the dashed histograms are the sum of the BDKRC 
simulated events and of the estimated background. The hatched areas represent the 
background contamination. 



DELPHI 






N 
•^ 



■§ 



3.5 r 



2.5 




1.5 - 



0.5 - 



Figure 5: The differential cross-section for the reaction 7*7* 
shown with error bars. The solid line shows the QED expectation. 
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Figure 6: Distributions of Ql (a), invariant mass Whad (b), charged multiplicity Ncharged 
(c), and Y calculated from the tagged particles' 4- momenta (d). The data are shown with 
error bars. The solid and dashed histograms correspond to the sum of 7*7* simulated 
events obtained with the PYTHIA and TWOGAM generators, respectively, and include 
the total background, shown hatched. 
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Figure 7: Reconstruction efficiency as a function of Y. The results of the calculations 
based on TWOGAM and PYTHIA generators are shown. 
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Figure 8: Differential cross-section for the reaction e"'"e^ ^ e+e^ + Ziarfrons . The dashed 
histogram corresponds to the average of the TWOGAM and PYTHIA predictions. The 
data are shown with error bars: the total error bars indicate the sum in quadrature of 
the statistical (inner error bars) and of the systematic uncertainties. 
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Figure 9: The differential cross-section for the reaction 7*7* -^ hadrons . The data 
are shown with error bars: the total error bars indicate the sum in quadrature of the 
statistical (inner error bars) and of the systematic uncertainties. The solid curve corre- 
sponds to the expectation of the quark-parton model (QPM, quark-box diagram, figure 
1). The two dotted lines represent the BFKL calculations in the leading order [12]. The 
next-to-leading order calculations [13] are shown by the two dashed curves in the mid- 
dle. The two curves for the BFKL calculations correspond to the Regge scale parameter 
changing between Q^ (upper line) and AQ"^ (lower one). The QPM contribution is added 
to both the LO and the NLO BFKL expectations. 



